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ON ROBUST SOLUTIONS TO UNCERTAIN LINEAR COMPLEMENTARITY 
PROBLEMS AND THEIR VARIANTS 

YUE XIE AND UDAY V. SHANBHAG* 


Abstract. A popular approach for addressing uncertainty in variational inequality problems is by solving the expected 
residual minimization (ERM) problem HCol. This avenue necessitates distributional information associated with the uncer¬ 
tainty and requires minimizing nonconvex expectation-valued functions. We consider a distinctly different approach in the 
context of uncertain linear complementarity problems with a view towards obtaining robust solutions. Specifically, we define 
a robust solution to a complementarity problem as one that minimizes the worst-case of the gap function. In what we believe 
is amongst the first efforts to comprehensively address such problems in a distribution-free environment, we show that under 
specified assumptions on the uncertainty sets, the robust solutions to uncertain monotone linear complementarity problem can 
be tractably obtained through the solution of a single convex program. We also define uncertainty sets that ensure that robust 
solutions to non-monotone generalizations can also be obtained by solving convex programs. More generally, robust counterparts 
of uncertain non-monotone LCPs are proven to be low-dimensional nonconvex quadratically constrained quadratic programs. 
We show that these problems may be globally resolved by customizing an existing branching scheme. We further extend the 
tractability results to include uncertain affine variational inequality problems defined over uncertain polyhedral sets as well as 
to hierarchical regimes captured by mathematical programs with uncertain complementarity constraints. Preliminary numerics 
on uncertain linear complementarity and traffic equilibrium problems suggest that the presented avenues hold promise. 


1. Introduction. The fields of robust control m and robust optimization j3] have grown immensely 
over the last two decades in an effort and are guided by the desire to provide solutions robust to parametric 
uncertainty. To provide a context for our discussion, we begin by defining a convex optimization problem 


min/(a;;M), (1.1) 

x^X 

where X C K”, u £ U C f ■. X x U ^ M. is a. convex function in x for every u £ U. The resulting 
collection of uncertain optimization problems is given by the following set: 

J ueu 

Given such a set of problems, one avenue for defining a robust solution to this collection of uncertain problems 
is given by the solution to the following worst case problem: 


min max f(x: u). 

xeX uGU 


( 1 . 2 ) 


Robust optimization has grown into an established field and there has been particular interest in deriving 
tractable robust counterparts to (1.2); in particular, can one formulate a single convex optimization problem 


whose solution lies in the set of solutions of ( |1.2[ ). Such questions have been investigated in linear, quadratic, 
and in more general convex regimes HIE] while more recent efforts have considered integer programming 
problems m- 

A particularly important class of problems that includes convex optimization problems is that of varia¬ 
tional inequality problems m- Recall that a variational inequality problem requires an x £ X such that 


{y - xY'F{x) >0, Vy e A, 


(1.3) 
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where F : X ^ M". Hereafter, this problem will be denoted by VI(X, F). Moreover when X is a cone, it is 
known that VI(X, F) is equivalent to the complementarity problem CP(X, F) (see [H]); The latter problem 
requires an x such that 

X^xFF{x)GX*, (1.4) 

where X* denotes the dual cone defined as X* = {y : y'^x > 0, x £ X} and y .L x implies that = 0 

for every i. Such problems have grown increasingly relevant in control and optimization theory and find 
application in the modeling of convex Nash games in communication networks Dun], traffic equilibrium 
problems m [19], and spatial equilibrium problems [22l |2T]. Naturally, in almost all of these settings, 
uncertainty represents a key concern. For instance, in Nash-Cournot games, the price function of the 
quantity being sold may have uncertain parameters while in traffic equilibrium problems, travel times are 
rarely known with certainty. Given such a challenge, one may articulate a parametrized VI(Jr, F(»; u)) that 
requires an x such that 

{y - x)^F{x;u) >0, Vj/G X, (1.5) 

where u denotes a belief regarding an uncertain parameter. The resulting collection of uncertain variational 
inequality problems is given by the following: 

{VI(X;F(.;w))}„^^. (1.6) 

In this paper, we consider the setting where F{x\u) = M{u)x + q{u) where M{u) G q{u) G K", and 

X = and the resulting collection of affine variational inequality problem is equivalent to collection of 
uncertain linear complementarity problems 


{LCP{M{u),q{u))}^^^ , 

(1.7) 

where hCP{M{u),q{u) is defined as follows: 


{M{u)x + q{u))^x = 0. 

(1.8) 

M(u)x + q{u) > 0, 


a; > 0. 



In subsequent sections, we extend these statements to uncertain affine variational inequality problems (AVI) 
over polyhedral sets and uncertain mathematical programs with complementarity constraints (MPCC). More 
general forms of F{x, u) will be considered in future research. Given that LGPs and MPGGs allow for 
capturing a broad class of application settings and the paucity of results for capturing uncertain forms 
of such problems, we believe that these questions have relevance and represent a necessary step before 
proceeding further. 

Now we briefly touch upon earlier efforts in addressing this class of problems. In particular, much of 
the prior work has focused on the minimization of the expected residual function (cf. [HI UHl HI] and the 
references therein). Given a random map F{x;^) where ^ H F : X x ^ K", and a probability 

space (H,F, P), based on the residual function f{x;^), the expected residual minimization problem (ERM) 
problem is defined as the following: 

minE[/(a:;^)]. (1.9) 

xGX 
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Such avenues have derived such solutions for both monotone as well as more general stochastic variational 
inequality problems but are complicated by several challenges: 

(i) First, such avenues necessitate the availability of a probability distribution P. 

(ii) Second, the expected residual minimization problem, given by ERM, leads to a possibly nonconvex 
stochastic optimization problem and much of the research has focused on providing estimators of 
stationary points to such problems. 

(iii) Third, this approach focuses on minimizing the average or expected residual and may be less capable 
of providing solutions that minimize worst-case residuals unless one employs risk-based variants. 

In the spirit of robust approaches employed for the resolution of a range of optimization and control-theoretic 
problems, we consider an avenue that requires an uncertainty set. An alternative of immense importance not 
considered here is the scenario-based approach [IlIH] and is left as the subject of future work. Furthermore, 
rather than minimizing the expected residual function, we consider the miminization of the worst-case 
residual over this uncertainty set. Specihcally, we make the following contributions: 

(a) First, in the context of stochastic linear complementarity problems with possibly asymmetric positive 
semidefinite matrices, we show that the robust counterpart is a single convex optimization problem 
under varying assumptions on the associated uncertainty set. 

(b) Second, we observe that somewhat surprisingly, robust solutions to non-monotone regimes can also 
be obtained tractably under some conditions, revealing that such problems are characterized by 
hidden convexity. More generally, robust solutions to uncertain non-monotone LCPs are shown to 
lead to low-dimensional nonconvex quadratically constrained quadratic programs. We customize a 
recently presented branching scheme to allow for obtaining global solutions to the resulting QCQPs. 

(c) Third, we extend our statements to two sets of generalizations: (i) Uncertain afhne variational 
inequality problems over uncertain polyhedral sets; and (ii) Mathematical programs with uncertain 
complementarity constraints. 

(d) These results are supported by preliminary numerics that provide a comparison across different 
solvers on a non-monotone uncertain TCP and a traffic equilibrium case study. In the instance 
of the latter, we observe qualitative distinctions between ERM solutions and the proposed robust 
solutions. 


The remainder of this paper is organized as follows. In Section we motivate our study through 
two applications and provide an instance of a monotone complementarity problem with arbitrarily high 
price of robustness. In Section we first define the robust counterpart and subsequenty discuss tractable 
reformulations of both uncertain monotone LCPs as well as a special case of uncertain non-monotone LCPs. 
Non-monotone generalizations and their global solutions form the core of Section while generalizations to 
variational inequality and hierarchical regimes are considered in Section]^ Preliminary numerics are provided 
in Section and we conclude with a short set of remarks in Section Finally, to improve readability, we 
list the key acronyms and their equation numbers: 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


VI or VI(A,F(.)): see 
CP or CP(A,E(.)): see Q. 

Uncertain LCP or uLCP{M{u),q{u)): see (1.8). 
ERM: see p.9|). 


MPCC: see (5.8). 
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2. Motivating examples and applications. In this section, we begin by providing an example of an 
uncertain linear complementarity problem where the lack of robustness has immense impact, and proceed 
to discuss two applications that motivate the study of uncertain linear complementarity problems. 


2.1. Robust solution vs. non-robust solution. We begin by considering the following simple 
uncertain LCP: 


“ V"'') 

\yj V® S{^,ri)J \yj 



V(u,^,77) 


( 2 . 1 ) 


where M = (^I - G S'(^, 77 ) = ^S'i+? 7 S' 2 ,x G K",?/ G G K!ji,g(w) = ue, u G lAu = [0,1]. 

Furthermore, Si = nin + and S 2 = + e„e^. e denotes the column of ones, e„ = (1,... and 

U(^ = {Q = ^ + r] >Q}. We begin by noting that a solution to the upper system 


0 < X _L Mx — Qx >0 

is uniquely defined by x* where x* = (/ + ee"’")qx > 0. The lower system requires solving the following 
equation: 


0 < y _L S{^,r])y + q{u) > 0. 


Since S{^, y) h 0 and q{u) > 0, it follows that y = 0 is a solution for all u £Uu and (5, y) G However, if 
^ = y = u = 0, then any nonnegative y is also a solution implying that there is a ray of solutions. Our focus 
lies in obtaining solutions that minimize the worst-case residual defined as follows: 


max H ((^ ° ]h] 



x^{Mx-qx)+ max iy'^ S{^,y)y + y'^ q{u)) 
«6Wu.({,OGWc 


= x'^{Mx-qx)+ niax Uy'^ Siy + yy'^ S 2 y + uy"^ e) 


= x^(Mx - qx) + maxly^S'iy, y'^S2y} + y^e. 


We use this setting to distinguish between a non-robust solution and a robust solution. 

A non-robust solution: Suppose the realization of m, y is such that m = 0,^ = 0,?7 = 0 and the resulting 
solution is given by (x*,yi) where yi > 0. Consequently, the worst-case residual given by x^(Mx — qx) -f 
maxjyJ’S'iyi, yf S' 2 yi} -I- yfe. could be arbitrarily high since yi is any nonnegative vector. In effect, a non- 
robust solution chosen under a single realization can have large worst-case residual. 

A robust solution: The robust solution of this problem is given by (x*,0) and achieves that worst-case 
residual equals to x*^(Mx* — qx)- 

Remark: This example shows that in contrast with a robust solution of an uncertain LCP, non-robust 
solutions may lead to arbitrarily poor worst-case residual. 

2.2. Applications. 

Uncertain traffic equilibrium problems. A static traffic equilibrium model [TS] HI] captures equi- 
libriating (or steady-state) flows in a traffic network in which a collection of selfish users attempt to minimize 
travel costs. Here, we present a path-based formulation in JV. JV denotes the network while A represents 
the associated set of edges. Further, let O and V denote the set of origin and destination nodes, respectively 
while the set of origin-destination (OD) pairs is given by W C O x P. Let Vw denote the set of paths con¬ 
necting each w GW and V = Liw^w'Pw Let hp denote the flow on path p gV while Cp(h-, u), the associated 
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(uncertain) travel cost on p, is a function of the entire vector of flows h = (hp) and the uncertainty u G U. 
Let diu(v,u) represent the uncertain travel demand between 0-D pair w and is a function of z; = (z’uj), the 
vector of minimum travel costs between any OD pair, and the uncertainty u G U. Based on Wardrop user 
equilibrium principle, users choose a minimum cost path between each 0-D pair: 


0 < hp ± Vw — Cp{h) >0, Vw e W, p G Vw 


( 2 . 2 ) 


Additionally, the travel costs are related to demand satisfaction through this problem: 


0 < -L ^ hp — dui{v) >0, Vw G W, 


(2.3) 


Static traffic user equilibrium problem is to solve a pair {h,u) satisfying (2.2) and (2.3), compactly stated 
as the following uncertain complementarity problem: 

^C{h;u)-B^v'^ 


0< (^1 _L I “ I >0, \/uGU 

\v I I Bh — d(v]u) j 


(2.4) 


where C{h]u) = {Cp{h]u) \ p G V), d{v,u) = (d„(z;;zt) | w G W) and B is the (OD pair, path)-incidence 
matrix {bwp)' 


h — 

^wp — 


1 if p GVw 
0 otherwise. 


This represents an uncertain collection of complementarity problems and we desire an equilibrium (/i, v) that 
is robust to uncertainty. 

Uncertain Nash-Cournot games. Nash-Cournot models for competitive behavior find application in 
a variety of settings, including the context of networked power markets [20] ■ We describe an instance of a 
single node A^—player Nash-Cournot game in which N players compete for a single good. Suppose player i’s 
uncertain linear cost function given by Ci(u)xi where Xi is her production level decision. Furthermore, the 
ith player’s capacity is denoted by capj. We assume that sales of the good are priced using an (uncertain) 
price function dependent on aggregate sales X and denoted by p{X;u) where u G U. We restrict our 
attention to settings where this price function is affine and defined as follows: p(A; u) = a{u) — b[u)X where 
o(m), b{u) > 0, A = uGU. The ith agent’s problem is given by the following: 


min {c;,{u)xi- p(X\u)xi) (Player(a;_i)) 

subject to Xi < capj. (A^) 


The sufficient equilibrium conditions of the Nash-Cournot game are given by the concatenation of the re¬ 
sulting optimality conditions: 


0 < T b{u){X -|- Xi) -I- -I- Ci{u) — a{u) > 0, 

0 < Ai T capj — Xi > 0, 

The resulting uncertain LCP is given by the following: 


Vi, 

Vi. 


0 < * ± A/(«) + ,(„) >0, V„ € U. where M{v.) 4 ' M , ,(„) 4 

\ -I 0 \ cap ! 


e denotes the column of ones, cap is the column of capacities, and 1 represents the identity matrix. 
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3. Uncertain LCPs with tractable robust counterparts. While optimization problems admit 
a natural objective function, such a function is not immediately available when considering variational 
inequality problems. However, one may define a residual/merit function associated with Y1{X^F). It may 
be recalled that a function g{x) is a residual/merit function for VI(X, F) if the following hold: 

(i) g{x) > 0 for a: G X; 

(ii) X € X solves VI(X, F) if and only if g{x) = 0 (see pi]b 

When X is a cone, the problem reduces to a complementarity problem and the residual/merit function 
associated with CP(X, F) is given by the following gap function, defined as follows: 


^gap 


{x,u) = 


F{x,u)'^x if F{x,u) G X* 


+00 


otherwise. 


(3.1) 


Throughout this section, we assume that the set X is the nonnegative orthant K", F{x,u) = M{u)x + q{u), 
where M(u) G and g(u) G K” for every u Gli. Furthermore, throughout this paper, we utilize the 

gap function as the residual function in developing tractable convex and relatively low-dimensional robust 
counterparts of uncertain LCPs with monotone and non-monotone maps. Specifically, in Section |3.1[ we 
consider settings where q[v) is an uncertain vector and M[u) is an uncertain positive semidefinite matrix with 
V GV and u G U. We provide robust counterparts in regimes where U and V are either distinct (unrelated) 


or related under varying assumptions on the uncertainty sets. In Section 3.2 we provide tractable robust 
counterparts to regimes where the M(u) is not necessarily positive semidefinite. 

3.1. Uncertain monotone LCPs. Much of the efforts in the resolution of uncertain variational in¬ 
equality problems has considered the minimization of the expected residual; instead, we pursue a strategy 
that has defined the field of robust optimization in that we consider the minimization of the worst-case gap 
function (residual) over a prescribed uncertainty set. While in its original form, such a problem is relatively 
challenging nonsmooth semi-infinite optimization problem. Yet, it can be shown that these problems are 
equivalent to tractable convex programs. By setting f{x,u) = 0gap(x, u), (0) can be recast as follows: 


min max Flx.u)^ x 
ueu 


subject to 


F{x, u) 

X 


G X* 
G X. 


yuGU, 


(3.2) 


Before proceeding, it is worth noting that the robust formulation attempts to find a solution that minimizes 
the maximal (worst-case) value taken by F{x, u)^x over the set of solutions that are feasible for every u GU. 


In fact, the following relationship holds between the optimization problem (3.2) and the original uncertain 
complementarity problem. 

Lemma 3.1. Consider the problem given by (|3.2[). Then x G X solves 


X ^ X 1. F(x, u) G X* for all u GU, 


if and only if x is a solution of (3.2) with optimal value zero. 


Unfortunately, it is unlikely that there exists an x that solves CP{X, F{»,u)) for every uGU\ instead, 
we focus on deriving tractable counterparts that produce global minimizers to ( |3.2[ ) which may be rewritten 
as follows: 


min max x'^{M(u)x + q{u)) 

a:>0 u£V( 


subject to 


min Mi,(u)x + qi{u) >0, VL 

u&A 


(3.3) 
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or 


min 

x>0,i 


subject to max {M{u)x + q{u)) <t, 

u&A 

min Mi,(u)x + qi{u) >0, Vi. 

u&A 


(3.4) 


We first consider the development of tractable counterparts of (3.3) or (3.4) when M{u) = M and q is 
uncertain. Subsequently, we consider the more general setting when both M and q are uncertain but are 
derived from unrelated uncertainty sets. Finally, we assume that both M and q are derived from related 
uncertainty sets. 

3.1.1. Uncertainty in q. In this subsection, q is subject to uncertainty and the problem can be reduced 
to obtaining a robust solution to an uncertain quadratic program (cf. [3]). We define q{u) as follows: 

L 

q{u) = qo + '^uiqi, u= {ui,U 2 , ■ ■. ,ul) GU, 

1=1 


(3.5) 


where U is a perturbation set yet to be specified. Consequently, (3.3) may be rewritten as follows: 

L 

min x'^{Mx + qQ) + maxy^ ui{x'^qi) 
u£U 




Vi 


subject to Mi,x + [go]i + min 7 ui[qi\i > 0 , 

U&J 

/ = 1 

X >0. 

We begin by considering three types of uncertainty sets: or U 2 -, where 

Ui^{u-. || m||i < 1 },U 2 = {u-. || m ||2 < 1}, and Uao = {u. || m||oo < 1} ■ 


(3.6) 


(3.7) 


The proof of our first tractability result is relatively straightforward and is inspired by Examples 1.3.2 and 
1.3.3 from [3]. 

Proposition 3.2 (TRC for uLCP{M,q{uj)). Consider the uncertain LCP given by (3.6) where M is 
a positive semidefinite matrix. Let U be defined as Ui,U 2 or Uao as specified by (3.7). Then (3.6) can be 
reformulated as a tractable program. 

Proof. 

(i) U := Uoo' We begin by noting that ||in||i = max Consequently, the term has 

lEIU^i 


maximal value Furthermore, we have that 

L 

'^ui[qi]i = 


mm 

U^lAcx. 


Z=1 


min u q = — max {—u q) = — 

UGUoo U£Uoo 


/ 91^ 

Vl/ 




z=i 


Consequently, (3.6) is equivalent to the following: 

min x'^(Mx + qo) + Jf;lL^lqfxl 
subject to + do — IqiI > 0 , 

X > 0 . 


(3.8) 
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Finally, by adding additional variables, (3.8) can be rewritten as a convex quadratic program (QP): 


min + go) + 

subject to ti > q[X 

Mx + qo- hi 

X 


> -ti, 
> 0 , 

> 0 . 


V/, 


(ii) 


U := Ui'. We proceed in a fashion similar to (i) and begin by recalling that 
leading to the folliwing simplification: 


max rFw = llw 
lhlli<i 


OO 5 


L 

maxY^ uix'^qi 

uGWi ^ 

1=1 




OO 


max \x^ qi\. 


Similarly, min rj^w = —||w||oo and 'Yld=iUi{qi)i has minimal value — max;g/i |(®)i|- Conse- 
quently, the constraints reduce to 


Mi,x + (qo)^ - max \[qi]^\ > 0, 

Vi 

Mi,x + {qo)^ - |(®)i| > 0, 

WiJ 

Mx + qo - hi > 0, 

yi. 


Similarly, the objective function can be stated as x'^(Mx + go) + max/g^^^ By adding a 

variable t, this problem may be reformulated as the following convex quadratic program: 


min x'^{Mx + qo) + t 



subject to t>q[x 

>-t, 

V/ 

Mx + qo- 1®! 

>0, 

V/ 

X 

> 0. 



(hi) U By leveraging Example 1.3.3 in 0, it is seen that 


T 

max n w = 
lhl|2<l 


T 

W W 

Iklb 


llw^lb- 


As a result, ^las maximal value while 


T 

mm Tj w = 
lhl|2<r- 


W^W 

Iklb 


-Ikll2, 


indicating that Yld=i ui{qi)i has minimal value — \JJ2iLihi)i ■ Then (3.6) may be rewritten as: 


min x'^ (Mx + qo) 
subject to Mx + qo — q 


> 0 , 

> 0 . 


X 















where qi = y J2iLi{Qi)i ■ By adding an extra variable, we obtain convex program with a quadratic 


objective and a conic quadratic inequality. 



min 

x'^{Mx + qo) + t 


subject to 


<t, 


Mx qo — q 

>0, 


X 

> 0. 


□ 

Next, we present a more result where the uncertainty set is captured by a more general convex set. Specifi¬ 
cally, U -—Uc where 

Uc= :3iy : Pu + Qi^+p€ K CW^} , (3.9) 


K is a cone, P and Q are given matrices, and p is a given vector. 

Proposition 3.3 (RC of uLCP(M, q(w)) when U := lAc)- Consider the uncertain LCP given by (3.6) 
where M is a positive semidefinite matrix. Let U := Uc, where Uc is defined as (3.9). Suppose one of the 
following holds: 

(i) K is a polyhedral cone; 

(ii) K is a convex cone and the following holds: 


3{u,v) such that Pu + Qv + p G int{K). 


(3.10) 


Then the robust counterpart of (3.6) is a finite dimensional convex conic program given by the following: 


min x"^M qx + qfjX 3- p^y 





subject to {P"'"y)i -\- x^qi 

= 0, 

yi = 

1,. 

..,L 

Q^y 

= 0, 




y 

e K*, 




Mi,x + [go]i -p'^Zi 

>0, 

Vi = 

1,- 

..,n 

[P,iVz, 

= fe]i, Vi=l,.. 

.,nj = 

1,. 

..,L 

Q^Zi 

= 0, 

Vz = 

1,- 

..,n 

Zt 

G K*, 

Vz = 

1,- 

..,n 

X 

> 0. 





The proof follows from Theorem 1.3.4 and Proposition 6.2.1 from [3] and is omitted. 

Remark: If K is chosen to be the nonnegative orthant, the uncertainty set is a polyhedron given <5 = 0. 
Both Ui and Uoo are included in this general case. If K is chosen to be the second-order cone, a special 
case of the perturbation set Uc is a ball. Under both circumstances, the problem is tractable. Notice that 


nonnegative orthants and Lorentz cones are self-dual. When K is chosen to be K" , (3.6) reduces to a convex 
quadratic program (QP). If K is chosen to be L”, (3.6) can be recast as a convex quadratically constrained 
quadratic program (QCQP). 


3.1.2. Uncertainty in both M and q under an independence assumption. Next, we consider the 
setting where both M and q are uncertain but the sources of uncertainty are independent. This is a somewhat 
more challenging problem and a direct application of the results from robust quadratic programming appears 
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difficult. 

Recall that the map F is said to be monotone over a set X if the following holds: 

{F{x)-F{y))'^{x-y)>0, \/x,y G X. 


Additionally, Mx + q is monotone over R" if and only if + M^) is positive semidefinite (cf. 115)1. 
Without loss of generality, we assume that M or XI(u) is sysmetric through this subsection, if not, we 
may always replace the matrices by their symmetrized counterparts. For the present, we assume that q is 
deterministic and reformulate (3.3) as follows: 


min t 

subject to x'^{M{u)x + q) < t, 
M{u)x + q > 0, 
a: > 0, 


Vu G U, 
Vu G U, 


where M{u) is defined as follows: 


L 

M(u) = Mq 

1=1 


Mo h 0, 

MiFO, l = l,...,L. 


(3.11) 


(3.12) 


Of course. Mi, I = 0,..., L have also been assumed to be symmetric. We now present a tractability result 
for nonnegative restrictions of and defined as follows: 


^oo = ■ Halloo < > 0} and Ui = {u : ||M|ji < l,u>0}. 


(3.13) 


Note that under the definitions of and U^, M{u) is always positive semidefinite. This implies that 3.11 
is convex for each u G U. 


Proposition 3.4 (TRC for uLCP(M(u),g) ior Consider the problem (3.11) where M{u) is 

defined by (3.12) and U is chosen either to be orUl, both of which are defined in (3.13). Then the 
uncertain TCP has a tractable robust counterpart, given by a QP and a QCQP, respectively. 

Proof. 

(a) U We first derive the robust counterpart of the following constraint: 


x^ Mqx F^^uix'^ Mix + x^q < t, yuGUf^. 


1=1 


But this can be equivalently stated as 


x^ Mqx + max 


uix"^ Mix 


L;=i 


+ x^q < t- 


We now evaluate the maximum in the right hand side: 


max 

u&AjL 


uix'^MiX 


. 1=1 


= > max > uiaMx^Mix, 0) = > x'^Mix, 

1=1 ‘ ^ ' 1=1 1=1 
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where the last equality is a consequence of applying the positive semidefiniteness of Mi for I = 
I,... ,L. Consequently, the robust counterpart of (3.11) can be stated as follows: 


min t 

subject to (Mq-\-Y^ i^-^Mi)x + q 


min 


X 


' Ui[Mi]t,x 


[Mo]i,x + qi 


<t, 

> 0, Vi 

> 0 . 


We may now simplify the second constraint as follows: 


where 


min 


■ L 

J^l 


L 


y min 


ui[Mi]i,x = 




[Mi]^,x, if [Mi]i,x < 0 

vi = < 

0, if [Mi]„x > 0. 


As a consequence, vi = — max(—0) for I = 1,...,L and the TRC may be rewritten as 
follows: 


mm 
subject to 


x'^{Mo+ Y^i^^Mi)x + x'^q <t, 
L 

-'^max{-[Mi]i,x,0) + [Mo]i,x + q^ > 0 , 

X >0. 


Vi 


1=1 


Through the addition of variables 

zi,... ,zl, we may eliminate 

the max. function. 

following quadratic program (QP): 




min 

x'^{Mo + J2t=i 



subject to 

Mqx + q- J2i=i zi 

>0, 



Mix + zi 

>0, 

yi, 


Zl 

>0, 

yi, 


X 

> 0. 



(b) U := U{-. In an analogous fashion, when U := Ul, we have the following sequence of equivalence 
statements for the quadratic constraint: 


x^MqX + max 


uix'^Mix 


. 1=1 


+ x^ q <t 



x’^MqX + 

max fmax(a;^M/x, 0)1 + x'^q < t 



ie{i,...,L} '■ ■' 


x'^MqX + 

max Ix'^Mix] + x'^q < t, 



/£{!,... ,Z/} 


x'^MqX + 

x'^Mix + x^q < t, 


1 = 1,...,L, 
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where the second equivalence statement follows from the positive semidefiniteness of Mi. The semi¬ 
infinite linear constraint can be reformulated as follows: 

r L 1 


min 


^Ui [Mi]i.x 
J^i 


+ [Mo]i,x + (7i > 0, 


Vi 




— max 


max{-[Mi]i,x, 0) + [Mo]i,x + qi>0, 


max{-[Mi]i,x,0) < [Mo]i,x + qt 
max{—Mix,0) < [Mo]x + q, 


Vi 

\/i,l 


Finally, by the addition of a variable z, we obtain the following QCQP: 


min 
subject to 


t 


x'^{Mo + Mi)x + x^q 

< 

t, yi 

Mqx + q — z 

> 

0, 

Mix + z 

> 

> 

o' 

z 

> 

0, 

X 

> 

0. 


□ 

Remark: Note that while we do not explicitly consider the case when q is also uncertain, this may be 
easily introduced when the uncertainty set that prescribes M{u) is unrelated to that producing q{u). On 
this occasion, we may address each term individually, as in the prior subsection. Next, we consider the 
possibility that M and q are derived from the same uncertainty sets. 


3.1.3. Uncertainty in M and q under a dependence assumption. Next, we extend the realm of 


applicability of the tractability result to accommodate uncertainty sets that are more general than (3.121. 
Specifically, we employ an uncertainty set that relies on computing the Cholesky Factorization of M, defined 
next as adopted in [^r 

UA={{M,q)\M = Vl^A, A = Ao + E 9 = 90 + E ’ 11^112 < 4 • (3-14) 

I i=i i=i ) 

We explore this construction since it allows for developing tractable robust counterparts when M{u) and 
q{u) are related. Consequently, (3.31 may be recast as follows: 


min t 

subject to x'^{Mx + q) <t, 
Mx + q > 0, 
X > 0. 


y{M,q)eUA, 

yiM,q)eUA, 


(3.15) 


For the sake of convenience, we write the first constraint as 

x^Mx + 2a;^ t < 0, V(M, q) e Ua 


(3.16) 


The tractability of constraint (3.16) follows from Theorem 2.3 in and is formalized without a proof. 
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Lemma 3.5. Consider the constraint (3.161 whereldA is defined by (3.14). Then the tractable counterpart 
of this constraint is given by (3.17); 


/ 

-PqX + t — t 

1 T 

1 T 

. , -^qfx 

(Aox) 


1 T 

-^qtx 

T 


{Aix) 


1 T 

-2qLX 


T 

{Alx) 

V 

AqX 

Aix, 

Alx 

I 


^ 0 . 


(3.17) 




However, it is more challenging to construct a robust counterpart of the constraint given by (3.18): 

Mx + q>0, y{M,q) gUa- 


(3.18) 


In fact, this is the key departure from the result provided in [^. For purposes of convenience and clarity, 
we rewrite Ua in terms of Aq, q^ and Ai,qi,l = 1,... ,L: 

CL L 

Ua = s {M, q) \ M = + ''^^{Af Ag + AqAi)S^i, + {Af Am + AmAi)f,if,m + Af Aiff , 


1=1 


Z<m 


/ = 1 


q = qo + 


1=1 ) 


(3.19) 


We may utilize [3l Lemma 14.3.7] in deriving the tractability of (3.18). 

Lemma 3.6. Consider the constraint (3.18) where Ua is defined by (3.19). Then the semi-infinite 
constraint has a tractable robust counterpart, which will be presented as (3.25). 


Proof. We begin by noticing that obtaining a feasible solution of (3.18) requires solving the following 
zth optimization problem for i = 1,... ,n: 

L L 

™ ^0 + Aq Ai]i,X + [qi]i) -G E +E 

1—1 l</<m<L 1 — 1 

subject to |j ^||2 < 1. (3.20) 


We may compactly rewrite (3.20) as follows: 

min 6i(x)^^ + ^^C'j(x)^ 
subject to ||^||2 < 1, 

where 6^ : M" —>• K^, : M" —>■ are all linear functions of x. We now define the following: 






,Mfix)=[ 




vC ' \kbiix) cfix) 


, and Z = 


' e 


II?I|2<1 


Then the QCQP (3.21) is equivalent to the following matrix optimization problem: 

min (^,Mi{x)) 

«62 


(3.21) 


(3.22) 


(3.23) 


where {A,B) = Tr{A'^B). Since the objective function is linear in we may extend the feasible region Z 
to its convex hull Z which is given Z = conv{2^}. By Lemma 14.3.7 from [3], we have that 


z = U = 


,w W 


e S 


L + l 


' 1 

,w W 


h 0,Tr(W) < 1 
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where represents the cone of symmetric positive semidefinite matrices. Using variable replacement, 

(3.23) is equivalent to the following semidefinite program: 


min {X,Mi{x)) 

pSi + l 

0 



< 1 , 



= 1 . 


(3.24) 


The feasible region of (3.24) obviously acquires a nonempty relative interior. Therefore its dual optimum 


can be obtained and thus we can reformulate the ith constraint in (3.18) as the following SDP constraint: 



(3.25) 


where ai{x) = [A'^AQ]i,x + [go]i- □ 

Our main result of this subsection can be seen to follow from Lemma 13.51 and Lemma l3. 6 


Theorem 3.7 (TRC for uLCP(M,( 7 ) for {M,q) € IAa)- Consider the uneertain LCP (3.15). Then this 
semi-infinite program has a tractable robust eounterpart given by the following SDP: 


min 

t,T,x,yi,y 2 


subject to T 





t 

+ Mo{x) >: 0, 

+ Mi{x) ^0, Vi = 1, 

OlxL^ 

yi,i + Via + > 0; Vi = 1, 

2/*.i <0, Vi = 1, 
x >0. 


,n 

,n 


where Mq : 

Remark: 


]^(L+r!,+ l)x(L+n+l)^ . 


R(.f-+l)x(-E'+l) Q,. 


are all linear functions of x. 


(i) When the uncertainty set ||^||2 < 1 is replaced by either ||^||i < 1 or ||C||oo < 1, Lemma 3.5 does not 
hold. What we may do instead is to enlarge the uncertainty set to get a tractable robust counterpart. 
In the case of ||^||oo < [31 Lemma 14.3.9] provides a semidefinite representable set that contains Z. 

On the other hand, we may enlarge ||Clli^lor||^||oo<lto their circumscribed spheres representing 


a scaling of U 2 allowing for the construction of tractable robust counterparts of (3.18). 
(ii) We note that 


claims a similar result (Theorem 3.2) as Theorem 3.7 However, there appears to 


be an issue in that the tractability of (3.18) is not proved and does not seem to follow directly. 


3.2. Tractable uncertain non-monotone LCPs. When the matrix ^{M{u)^ + M(m)) is not neces¬ 
sarily positive semidefinite for every u G U, M{u)x -I- q{u) is no longer monotone for every u. Consequently, 


the problem (3.11) no longer has convex constraints for every realization of u. As we proceed to show, we 


may still obtain a tractable robust counterpart under a suitably defined uncertainty set on M{u) with the 
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caveat that M{u) is unrelated to q{u). We begin by defining the uncertainty set for M{u). 


M{u) = Mo + 


1=1 


-(Mo + Mo^)^0, 


(3.26) 


Without loss of generality, we assume that M{u), Mq and M; are symmetric for Z = 1,..., L; if not, we may 
always replace the matrices by their symmetrized counterparts. Thus far, the definition of M(u) is the same 
as (3.12). However, we will takelA osUx^Ui^Uao- This approach leads to indefinite matrix thus nonmonotone 
affine map for some scenarios. The tractability of the robust counterpart of the uncertain nonmonotone TCP 
is proved next. 


Proposition 3.8 (TRC for non-monotone uLCP(M(it), g)). Consider the problem (3.11). Suppose M 
is defined by (3.26) andU is either Uaa Mi orlA 2 - Then this problem admits a tractable robust counterpart. 
Proof. 

(a) U := {u : ||u||oo < 1}: We begin by determining the robust counterpart of the following constraint: 

L 


x"'"M qX+ ^^uix^Mix + x^'^q < t, yuGlt. 
1=1 

This may be equivalently stated as 

x’^MqX + max 


(3.27) 


U^UcK 


uix^MiX 


. 1=1 


+ x^q < t. 


By noting that the summation can be written from I = 1,..., K, through the application of 
max||„||^<i rj'^u = ||??||i, it follows that 

L 

x’^Mqx + \x'^Mix\ + < t. 

1=1 


Since Mi P 0 for every I, \x'^Mix\ = x'^Mix. Consequently, the robust counterpart of (3.27) can be 
stated as the convex constraint: 

/ L \ 


Mo + '^Mi\ x + x^q < t. 


1=1 


Similarly, the constraint M{u)x -I- g > 0, \lu^U can be reformulated as follows: 
M{u)x + q>0, yuGU 




Mqx + uiMix -I- g > 0, Vu G 77 Mqx + min 


1=1 


Mqx — max 

U&Aoo 


y^ Ui [-Mix] 
. 1=1 

Mqx — \Mix\ -I- g > 0. 


L 


y^^uiMi3 


.1=1 


+ q>0, 


+ q> 0 Mqx — I — Mix\ -f g > 0, 


1=1 


1=1 
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Through the addition of variables, zi,..., zl, the resulting robust counterpart can then be stated 
as the following convex QP: 


(3.28) 


min x'^{MoMi)x + x'^q 




subject to MqX + q — Yld-i 

>0, 



zi > Mix 

> -zi, 

1 = 1,. 

..,L 

Zl,X 

> 0. 

1 = 1,.. 

.,L 


(b) U :=Ui As in (a), we begin by determining the robust counterpart of (3.27|: 

■ L 

uix^ Mix 


X MqX + max 

u^Ui 




+ < t. 


By noting that the summation can be written from I = through the application of 

max||i;|p<i = ||f?||oo, it follows that 

x^Mqx + max Ix'^Mixl + x'^q < t. 

But x'^Mix > 0 for all x and I = 1,..., L implying that this constraint can be rewritten as follows: 

x'^Mqx + max {x^Mix} + x^q < t. 

The max. function can be eliminated by replacing each constraint by a finite collection: 

x'^M qx + x'^Mix + x'^q <t, I = I,..., L 

Similarly, the semi-infinite constraint M{u)x -I- <? > 0, Wu GU can be reformulated as follows: 


M{u)x + q>0, Mqx+ ''^^uiMix + q>Q, \/u&U 


Mqx + min 
u&A\ 






1=1 


q>0 ^ Mqx — max 

uGUi 


+ q>0, 


Mqx — max [I — Mix\] + q > 0 
Consequently, the TRC is given by the following: 


ui [-Mix] 

. 1=1 

Mqx — z + q > 0, 
z>Mix>—z, V/ = 1,...,L^ 


min t 


subject to x^{Mo + Mi)x + x^q <t, \/l 

Mqx + q — z > 0 , 

z > Mix > —z, \/l 

x,z >0. 


(c) U :—U 2 - We first consider constraint (3.27) which can be equivalently stated as follows: 


x^ MqX + max 
uGU2 


y^ uix'^Mix 


. 1=1 


+ x^q < t 


X^ MqX + 


\ -|- x'^q < t. 

\ 1=1 
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Similarly, the constraint M{u)x + 9 > 0 for every u GU = U 2 can be reformulated as follows: 


M{u)x + g > 0, 


'iuGU 


^ MqX + min 
U&A 2 


MqX — max 
U&A 2 


^ uiMiX 


.1^1 
■ L 


+ g > 0, 


^ui [-Mix] 


. 1^1 


+ ^>0, 


■O' 


, ^ + gi > 0, Vb 

\ 1=1 


Consequently, the robust counterpart of (3.11) can be stated as: 


subject to x'^M qX + \JMix)"^ + < t 

[Mq]„x - y^Ef=i [|[^i]«a;|]^ + q^ > o, 

X >0. 

By examiniirg the second derivative of /(x) defined as 


/(^) = 


L 


yi 


(3.29) 


. ^(x^M/x)2, 

\ 1=1 

it can be concluded that / is a convex function. This result indicates that the the left hand side 


of the first constraint in (3.29) is a convex function, implying that the resulting feasible region 


is convex. The n remaining inequalities are in the form of second-order cone constraints and are 


therefore tractable convex corrstraints. It follows that (3.29) is a convex program. 


□ 

To get a geometric understanding of the prior proposition, we consider the following example. 
Example: Consider the case when M(u) and g are defined as follows: 


M{u) = u 


1 O' 

,0 2 , 


, g = I 1 , and U = {u \ —1 < u < 1}. 


It can be observed that the constraint: x'^{M{u)x -I- g) < t is nonconvex when u < 0. Note that this 
constraint can be rewritten as u{x\ -\- 2x1) -I- 2xi -I- 2 x 2 < t and we defined Ru as follows: 

Ru — {{xi,X2,t) I u{x1 + 2x1) + + 22^2 < t} ■ 


Then i?_i denotes the region above the surface shown in Fig. 3.1 labeled u = —1, clearly a nonconvex 


set. Likewise, the feasible regions Ri,Rq represent the regions above the surfaces presented in Fig. [3T 
labeled u = 1 and m = 0, respectively. Though the set R-i is nonconvex and appears to make the program 
challenging to solve, a better understanding emerges when we consider the intersection of R^ over u, as given 

by 




n R. 

-l<u<l 
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It can be seen that i?i = R. The figure on the left in Fig. 3.1 hints as to why this holds. The three surfaces 
intersect at a single point, namely (0, 0) and the surface with greater index u stays above that with the lower 
index. This implies that C i?Q C R_i. Actually, i?„ is monotone in u in that if ui < M 2 , then ^ Ru^- 
When considering such constraints in higher dimensions, similar behavior emerges. Finally, there have been 
prior observations regarding the presence of hidden convexity in nonconvex programs (cf. [3]). 




1 

0 

-1 


Fig. 3.1: Hidden convexity in two and three dimensions 


3.3. A general tractability result. In the prior subsections, we have provided a tractability statement 
when M(u) and qiu) are defined as per 

L L 

M{u) = Mo + ^w/M/ and q{u) = -^y^mqi. 

1^1 1^1 

A natural question is whether a more general tractability statement is available. We address precisely such 
a question by assuming that (M, q) S W, where W represents a more general uncertainty set. Note that in 
the settings considered in the prior snbsection, W takes the form given by the following: 

{ L L 

{M{u),q{u)) : M{u) = Mq + '^ uiMi and q{u) = <?o + ^ ^ 

1=1 1=1 

We address the tractability question by considering the related separation problem. This requires a crucial 
result that relates the tractability of the separation and the optimization problems, both of which are defined 
next over a compact convex set. 

(OPT) Optimization problem: Given a vector c £ Q" and a number e > 0, find a vector y G Q" such 
that y is an e-feasible and e-optimal solution of the problem: 

T 

max c X. 

If d{y,Z) denotes the Euclidean distance of y from Z, then y satisfies the following: 

d{y,Z) < e (e-feasibility) and c^x < c’^y -I- e, for all x G Z (e-optimality) 
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(SEP) Separation problem; Given a vector y G Q and a parameter e > 0, one of the following may be 
concluded: 

(i) assert that d{y,Z) < e; 

(ii) find a vector c S Q” such that ||c|| > 1 and for every x G Z, c^x < c^y + £■ 

Then Grotschel et al. m showed that these two problems are equivalent through the following result. 
Lemma 3.9 (Grotschel et al. |18|L Suppose Z is bounded with a nonempty relative interior. Then 
there is a polynomial algorithn^for solving the separation problem for Z if and only if there is a polynomial 
algorithm for solving the optimization problem for Z. 

In the context of obtaining a robust solution to an LGP, the optimization problem of interest is given 
by the following: 

min cf" z 

{ x^{Mx + q) <t, V(M, q) GW 
z={x;t): Mx + q>0, \l{M,q) GW 

x>0 

where d = (0; 1). 

Theorem 3.10 (Tractability of RC to uncertain LCPs with general uncertainty sets). Suppose 
there exists a polynomial algorithm to solve the separation problem over the compact convex uncertainty set 
W and the following assumptions hold: 

(a) M h 0,V(M,g) e W. 

(b) If denotes the ith row of M, then for i = I,... ,n, the following holds: 

||M,,.|i 2 > AeQ++, y{M,q)GW. 



(c) Both W and Z are bounded sets, each having a nonempty relative interior. 


Then the problem (3.30) may be solved by a polynomial algorithm. 


Proof. It suffices to prove is that separation problem over the set Z can be solved in polynomial time 
since the cost vector d in (3.30) has rational entries; specifically, d = (0,... ,0; 1)^. 


First notice that the assumption (a) implies that Z is convex since it is defined as the intersection 
of an infinite collection of convex sets. Therefore, the separation problem over Z is defined in a general 
form by (SEP). We proceed to show that this separation problem may be solved in polynomial time. It 
suffices to show that either (i) or (ii) in (SEP) can be shown to hold in polynomial time for a given point 
z = {x-,t) G and a positive parameter e. 

Gonsider the feasibility of z with respect to Z. Given a vector z = {x;t) where x is nonnegativ^ it 
suffices to examine whether z is feasible with respect to the following: 


x^{Mx + q)<t, y{M,q) gW = max x'’"{Mx + q)<t (3.31) 

{M,q)GW 

For f = 1,..., n, (Mi,x + qi)>0, y(M,q) GW = min (Mi,x + qi)>0, (3.32) 

’ (M,q)€:W 


^The algorithms are polynomial with input length n + log(i). 

^Note that if x has any negative component, a separating hyperplane can be constructed with relative ease and we ignore 
this possibility without any loss in generality. 
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In effect, feasibility can be ascertained if the following problems may be solved in polynomial time: 


(1). < max x^Mx-\-q^x> and (2i). < min Mi,x + qi>, i = 

{{M,q)ew j l(M,g)ew ’ j 

Given an (x; i) € problem (1) is linear in (vec(-M); q) and can be recast as a problem of maximizing 

a linear objective over a convex set. Furthermore, an ei-optimal solution of (1) will be denoted by (Mo,( 7 o) 
where ei will be subsequently defined. Similarly, problem (2i) is a convex program with an objective that is 
linear in {Mf^]qi) and thereby linear in {M;q). An e 2 -optimal solution of (2i) will be denoted by 
where £2 will be subsequently defined and the ith row of M® and g® will be collectively denoted by {ai,bi). 
Since the coefficients of both cost functions are rational by assumption, (1) and (2i) for i = 1, ... ,n can be 
solved in polynomial time upto any precision ei > 0 and £2 > 0, respectively and provide rational solutions. 

By the boundedness of Z, we may assume that ||y|| < fx for all y € Z. Prior to proceeding, we quantify 
the precision levels £1 and £ 2 . Given an £, suppose £1 and £2 are defined as follows: 

£1 A min{£/(3||a;|p + 3||x|| + 1), £/(3(^ + ||a:||)(2||x|| + 1))} and £2 = min{A£/( 6 /x), A£/ 6 , A/2}. 

Suppose yVi = {(Mi_,,gi) : {M,q) € W}. Since d{{Mo,qo),}V) < £ 1 , d{{ai,bi),Wi) < £ 2 , and W is com¬ 
pact, we may find (Mi, gi) € W and {di, bi) G Wi such that d{{Mo, go), (Mi, gi)) < £1 and d{{ai, bi), {di, bi)) < 
£ 2 . Hence, 


\x'^{Mox + go) - x'^iMix + gi)| = |a;^((Mo - Mi)a; -f (go - gi))| 

< ||a;||||(Mo-Mi)a;-b(go-gi)|| 

< ||a;||(||Mo - Mi||||x|| -b ||go - gi||) 

< ||a;||(||x||£i -b£i) 

<(||xf+ ||:r||)£/(3||a:f+ 3 ||x|| + 1)<£/3. (3.33) 


Furthermore, we have the following: 

VyeZ, we have \{{Mo + Mq )x + qo)'^{y - x) - {{Mi + Mf )x + qi)’^{y - x)\ 

< \\y - a;||||(Mo - Mi + - MI)x -b go - gi|| 

< ||y- a:||(2||Mo - Mi||||x|| -b ||go - gi||) 

< ||y-a:||(2||a;|| -b l)£i 

<{d+ l|a;||)(2||a;|| -b l)£i < e/3. (3.34) 

Similarly, we may bound the difference between {ai,bi) and {di,bi) for all v. 

for f = 1,..., n, we have |6i/A — 5i/A| < £2/A < e/6 (3.35) 

and for z = 1,... ,n and for all y £ Z, we have |(ai/A — di/Xj'^y] < ||ai — di||||y||/A 

< Iki - di\\y/X 

< Me2/A 

< e/ 6 . (3.36) 


20 



In addition, since problem (2i) is solved to a precision of £ 2 , we have that the following holds: 


I llozll - lla^ll I < £2 < A/ 2 . 


Then we may conclude the following: 


=>||a*||/A > ||a,||/A - 1/2 > 1 - 1/2 = 1/2, (3.37) 

where the last inequality is a consequence of the assumption ||ai|| > A. 

The final part of the proof is reliant on considering three possibilties that emerge based on obtaining 
the solutions to problems (1) and ( 2 i) for i = 1 ,..., n. 

(1) Suppose MqX + q^x < t — ei and UiX + > £2 for i = 1,..., n. Since (Mg, go) is an £i-optimal 

solution, we have that 

max (x'^Mx + q^x) < x'^MqX + gJa; + ei < t — £1 + £1 < t, 

(M,q)eW 

where the second inequality follows by our assumption on (Mg, go). Similarly, we have that {Mi^,,qi) 
is an £ 2 -optimal solution implying that 


min {Mi ,x + qi) > atx + - £2 > £2 - £2 = 0 , 

(M,q)GW 

where the second inequality again follows by assumption. It can be concluded that {x; t) G Z and 
(i) of (SEP) has been concluded in polynomial time. 

(2) Suppose x'^Mqx + q^ x > t — ei and UiX + > £2 for i = 1,... ,n. Then we may construct a vector 

c defined as c = [(Mo + Mq)x + go; —1]. Consequently, ||c|| > 1. Furthermore, since Mg and go have 
rational entries by assumption, it follows that c G Q”+^. Furthermore, V(j/;t) feasible, we have the 
following sequence of inequalities: 

r > y^{Miy + gi) > x'^{Mix + gi) + [(Mi + MI)x + qi]’^(y - x), 


where the second inequality follows from the convexity of y^{Miy + gi). It follows that 


x'^{Mix + gi) + [(Ml + mI)x + qiY'iy - x) 

> x^(Mox + go) - e/3 + [(Mi + Mf)x + gi]^(y - x) 
>t-ei-e/3+ [(Mg + Mq )x + go]’" ( 2 / - x) - e/3 
= t-e + [(Mo + Mg)x + qo]^(y - x). 


where the first inequality follows from (3.331, the second inequality follows by assumption on 
x’^Mqx + q^x and by invoking (3.34). Therefore we have that 


c^z = [(Mo + Mq)x + go; -l]’’(x; t) + e 

> [(Mo + M^)x + go; -l]’’(y; r), V(y; r) e Z. 


Consequently, we have concluded (ii) of (SEP) in polynomial time. 
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(3) Suppose there exists an i G n} such that ajx + h < € 2 - Then let c be defined as c = 


— (2ai/A;0). Then from (3.37), ||c|| = ||2ai||/A > 1 and c G Q"'+^ since G Q", A G Q. Thus 


V( 2 /;t) G Z,c^(y;T) = -2aJy/X 

< —2aJy/\ + 2e/6 

< 26,/A + e/3 

< 26,/A+ 2e/6 +e/3 

"G — 2(1^Xj\ -\- 2e2/A -\~ 2e/3 

< —2ajxl\ + e = <P"{x] t) + e, 


where the first inequality follows from (3.36), the second inequality follows from the feasibility of 


(hi, 6i), the third inequality is a consequence of (3.35), the fourth inequality is a consequence of the 
assumption that af x + bi < 62 , and the last inequality is a result of invoking the definition of e 2 . 
Therefore, we have concluded (ii) of (SEP) in polynomial time. 

Based on (1), (2), and (3), we note that the separation problem (SEP) can be solved in polynomial time. 
Consequently, by Lemma |3.9| there is a polynomial time algorithm for solving the optimization problem 
over Z. □ 

Remark: Through the above proposition, we establish a connection between the tractability of the 
uncertainty set and the tractability of the robust counterpart of uncertain LCP, further generalizing our 
findings from the prior subsections. 

4. General uncertain non-monotone LCPs. In this section, we consider non-monotone uncertain 
LCPs in more general settings where tractable robust counterparts are unavailable. Instead, we examine 
when such problems result in finite dimensional nonconvex programs. We assume that M{u) takes a form 
given by M{u) = Mq uiMi and q{u) = qo + where u G U, (Mg, go) denote the nominal 

values on M and g, and {Mi,qi)f^-^ are the basic shifts, while U represents a given perturbation set. This 
model of utilizing nominal values and shifts has been relatively standard in the field of robust optimization 
(cf. [3] and the references therein). In Section 4.1 we see that in the more general setting, the RC is 


already intractable to resolve when the perturbation set is of dimension one, demonstrating the hardness of 
the problem. In Section [4.2[ we discuss two avenues via which we may obtain nonconvex robust counterparts. 
While stationary points of such problems can be obtained by nonlinear programming solvers, global solutions 


require branching-based schemes. In Section 4.3 inspired by recent research by Fampa et al. m, we present 
a technique for obtaining global solutions of a nonconvex quadratically constrained quadratic program. 

4.1. NP-hardness of the problem. Throughout this section, we define M(u) as follows: 

L L 

M{u) = Mo + ^ ^^mquu € U. (4.1) 






without any assumption on Mg or Mj, / = 1, ..., L. Next, we prove that even the RC of a simple LCP(M(u), q{u)), 
where M{u) = Mg -|- uMi,q{u) = go -I- ugi, u GU = [Q,l] is NP-hard. 

Lemma 4.1. Consider an uncertain linear complementarity problem LCP{M(u),q[u)) where M(u) = 

Mg -|- uMi, q{u) = qo + uqi, u G U = [0,1] and Mg and Mi lie in K"^". Then the robust counterpart of this 
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problem is NP-hard where the robust counterpart is defined as follows: 

min max {M{u)x + q{u) 

u&A 

subject to M{u)x + q{u) >0, Mu^U (RClcp) 

X >0. 

Proof First we write the RC, given by (RCiCp), as follows: 

min max x’^ {M{u)x + q{u) 

u&A 

subject to M{u)x + q(u) >0, 'iu^U (RClcp) 

X >0. 

Let Mo = 0 and qg = 0, then 

max{a;^(M(u)a; + q{u))} = max {x"^{uMix + uqi)} = max{0, x'^{Mix + gi)}. 

ueU uG[0,l] 

Furthermore, we have that 


M{u)x + q{u) > 0,Vu S [0,1] Mix + > 0. 


As a consequence, max{0, x'^{Mix + gi)} = x'^{Mix + qi) for any nonnegative vector x and the (RClcp) is 
equivalent to the following: 


min 
subject to 


x’^{Mix + q) 

Ad\X ~\~ q\ ^ 0 ? 
X > 0, 


which is essentially the problem LCP(Mi, qi). But LCP(Mi, qi) is NP-complete since the equality-constrained 
0-1 knapsack problem can be reduced to LCP(M, q) in polynomial time [12]. Consequently, obtaining a ro¬ 
bust solution to an uncertain LCP with general matrices is NP-complete. □ 

Remark: If we assume that the feasible region {x \ Mix -I- qi > 0, a; > 0} is bounded, the LCP(Mi, qi) 
is still NP-complete. The reason is that equality-constrained 0-1 knapsack problem can be reduced to an 
LCP(M, q) such that {a; | Mx -I- q > 0, a; > 0} is bounded. Please refer to [12] for details. 


4.2. Nonconvex robust counterparts. In Section [3l2] we showed that the RC of uncertain nonmono¬ 
tone LCP may be tractable under some assumptions. However, in general, this is not the case, particularly 
when M(u) and q{u) are defined on the same (or related) uncertainty sets. In this circumstance, we can 
show that the RC may still be reformulated as a finite dimensional mathematical program. 


Proposition 4.2 (For non-monotone uLCP(M(u), q(u))). Suppose M(u) and q{u) are defined as (4.1) 
where U is defined as hli,hl 2 : or Uoo, (fdoo is defined in ( [3.7 ] ) ). Then ( [3.4] ) may he written as a finite 
dimensional mathematical program. 

Proof. We introduce an artificial variable w into a finite set of nonconvex quadratic equality constraints: 


wi = x'^Mix + qf X, Z = 1 ,..., L, 


Then the quadratic constraint in (3.4) can be written as follows: 


x^Mqx q^x + u^w < t, 


'iu&U 


(4.2) 
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(4.3) 




This semi-infinite constraint can be equivalently written as follows: 


Mqx + Qq X + max w <t. 
u^U 


(4.4) 


But, maxu^ui is given by ||w||i, ||rc||oo or ||rc ||2 if is given by UooMi or U 2 , respectively. Consequently, 

u&A 

(4.3) may always be recast in a tractable fashion. Similarly, consider the second constraint in (3.4): 

L 

Mqx + qQ+'’^ui{MiX + qi) >0, yu&U. (4.5) 


Based on Prop. 3.8 the constraint (4.51 can also be reformulated in a tractable fashion when U = U\,U 2 
or Uaa- We demonstrate this in detail when U = Uao and omit the proofs when U =U\ and U 2 - Note that 

(4.6) 


when U — Uoo-, (4.4) can be rewritten as 


X MqX + qQX+ llmlli < t. 


Furthermore (4.5) is equivalent to 


(Mox), + iqo)^ - ll^^lli > 0, {Z^)i = {Mix + qi),, W = 1,..., L, Vt = 1,..., u. 


Through the addition of variables, ti, ... zi,..., zl, the resulting robust counterpart when U = Uoo can 
then be stated as the following optimization problem: 


min 

x'>Q,t 

t 



subject to 

—Ti < x'^Mix -j- qf x 

<Tl, 1 = 1,. 

..,L 


x'^Mqx -i-q’^x-{- n 

<t, 



MqX -\-qo- 

>0, 



zi > Mix -1- qi 

>-Zl, 1 = 1,.. 

..,L. 


This QCQP is nonconvex when Mi, I = 1,... ,L are indefinite. □ 

We now provide a corollary of this result when U ^Uc- 

Corollary 4.3 (RC for non-monotone uLCP(M(u),q(u)) where U = Uc)- Suppose M{u) is defined 
as (4.1) and q{u) = go + Sfci where u£U andU is given byUc, defined as (3.9). Then (3.3) can be 
represented as a finite-dimensional mathematical program: 


min 

ai>0 

x’^Mqx q^x p^y 




subject to 

{P^y)i + x'^Mix qf X 

= 0, 

yi = i,.. 

.,L, 


Q^y 

= 0, 




y 

e K*, 




[Mo]i,x -1- [qo]i -p^Zi 

>0, 

yi = i,. 

..,n 


\P,lV Zi 

= [Mi]„x-i-[qi]i, Vi=l,.. 

.,n,l = l,. 

..,L 


Zi 

= 0, 

Vi = 1,. 

..,n 


Zi 

e K*, 

Vi = 1,. 

..,n. 


Remark: While stationary points of such problems may be computed through standard nonlinear pro¬ 
gramming schemes such as globalized sequential quadratic programming or interior point methods [26) . 
our interest lies in obtaining global solutions of such problems. In the next subsection, we review several 
approaches for obtaining global solutions to QCQPs. 
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4.3. A branching scheme for resolving nonconvex QCQPs. Before presenting our scheme, we 
provide a brief review of global optimization schemes for resolving indefinite quadratic programs and their 
quadratically constrained generalizations. Such a class of problems has seen significant study [a [23]. In la, 
the authors combine reformulation-linearization-technique (RLT) with an SDP relaxation to tackle QCQP. In 
[23| . a general framework is buit for solving such problems. While branching schemes come in varied forms, 
Burer and Vandenbussche [B] employ SDP relaxations for addressing indefinite quadratic programming. 
We consider a spatial branch-and-bound approach inspired by Fampa et al. [T6| developed for nonconvex 
quadratic programs. This approach uses secant inequalities for deriving a relaxation of the quadratic objec¬ 
tive. We extend this approach to quadratically constrained variants. We emphasize that the focus of this 
paper lies in extending standard robust optimization techniques to allow for accommodating uncertain linear 
complementarity problems. While a comprehensive study of branching schemes is beyond the scope of the 
current paper, we show that at least one of the approaches can be readily adapted to this context. 

We begin by noting that the QCQP can be recast as an optimization problem with a linear objective 


and quadratic constraints, some of which may be nonconvex. We continue using the defnition (4.1) and 
illustrate the scheme for the case when U := = {u \ ||u||oo < > 0} and qualify the relaxations and 

the bounds by using the superscript oo. Suppose u G W^. From Prop 14.21 the optimization problem given 


by (3.4) may be reformulated as follows: 


mm 


subject to x'^Mqx + QqX + X]i=i xaax{x'^MiX + qjx, 0} < t, 

{MQX + qo)k + YA^iX[m\{{MiX + qi)k,0} >0, VA: = I,...,n. 


(4.7) 


While the second row of the constraint part in (4.7) can be immediately written as: 

L 

Mox + qo-^z^>0, MiX + qi + Zi>0, z^>0, Vi = I,...,L, 

the chief concern lies in the first constraint which can be decomposed into and 2L + 1 constraints: 

L 


X 


Mqx QqX + ^ Ti < t, and > 0, x"^Mix + x < r*, z = 1,..., L. 


(4.8) 


i=l 


Constructing a relaxation: Akin to the approach employed in [T6|, we use the eigenvalue decomposition 
of Mi, defined as 

Ji Ki 

where > 0,Vj,Vi = I,...,L. Let M+ = = 0, ...,L. By dehning as 


ViJ ~ quadratic inequalities in (4.8) may be rewritten as 

Jo L 


i=i 


“ X! + ^’^^0 x + qoX + ^n<t, 

i=l 

Ji 

^ + Ji ^ i = I,...,L. 


(4.9) 

(4.10) 
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Suppose lij < Uij < Uij for j = and i = 0,...,L. Then we may use a secant inequality for 

providing a relaxation to the (4.9) in the form of the following: 


Jo L 

- X! ((yoj “ (“Od + ^Od) + ^0,j) + x + q^x + '^Ti^t, 
i=i i=l 


Ji 

~ X! ('“dj + ^ij) + ^ij) + x^M^x + qJx<Ti, i = 1,..., L. 

i=i 


(4.11) 

(4.12) 


When Z// := the resulting relaxed problem is denoted by {Pod{1,u)) and is defined as follows: 


min t 


{ivoj — h,j) (“o,i + h,j) + ^oj) + x'^Mq X + qQX + A 

n 

IV lA 
0 

z = 1,.. 

.,L, 

— 'Ylij=i {iViJ ~ ^i,j) A A h,j) A X Afj X + q^ X 

< Ti, 

z = 1, .. 

.,L, 

VI 

— ’ 



Vi ,3 - '/Kj’^IjX 

= 0, 

yj,h 


MqX + qo- X)f=i 2:i 

MiX + qi + Zi 

>0, 

>0, 

Vz = 1, . 

,..,L 

Zi 

>0, 

Vz = 1,. 

,..,L 

X 

> 0. 




{Poo{l,u)) 


Obtaining upper and lower bounds for yij: Crucial to this scheme is the need for obtaining upper and 
lower bounds on yij given by Uj and utj, respectively. Consider the set X^o 


Y ^ 


3zi G M”, Mox + go - ^-Zi > 0, MtX + qi + Zi>0, Zi > 0, Vi = l,...,L, a; > 0 


The lower and upper bound for yij can then be obtained by solving the following set of linear programs: 


• / - T 

min / max y 

subject to X G X^. 




Note that we assume these LPs are bounded and when the uncertainty set is either Ui or U 2 , the relaxation 
and the upper/lower bounds have to be derived in an analogous fashion. A formal outline of the branching 
scheme is provided in Algorithm]^ 
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Algorithm 1 Spatial branch and bound 


1; Init.i: k := 1; terminate := 0; choose e > 0; 

2: Init.ii: Mi = - + Y.fli for i = 0,..., L; 

3; Init.iii: For all i,j compute Ifj and ; 

4 ; Init.iv: Let Pi := Pacil°°,yl, G argminPi; 

5; Init.v: Assign bounds: glb_lb := t*; glb_ub := xl'^M qx^ + q^x^ + J2i=i ^SLx{xf MiX^ 
6: Init.vi: Update list: list := {(Pi,glb_lb,glb_ub)}. 

7: while terminate == 0 do 

8: Branching index: For Pk, choose index pair (z, j) := argmax(Mij — kj)-, (l>ij = j 


9: Update bounds: 


+ qlx*i,0}; 

+ j)/2; 





(bj) = (bj) 

otherwise. 


and lij 


(bj) = (bj) 
lij, otherwise. 


10 : Construct leaves: Pj. Pao{u,l); P^ := Poo{u,l). 

11: Upper and lower bounds for Pj.: {xl,Tl,yl, G argmin(P;J); 


L 

lb' ub' := {x^Y'MqxIp q^xi + ^ max{(xl)'^M,a;t + gf 0} 

i=l 

12: Upper and lower bounds for P^: {x^,T^,y^,z^,t^) G argmin(P^); 

L 

lb“ := ub“ := {x'^Y'M qxI Pq^x^ + ^ max{(x“)'^Mia;“ +<7^2;“}; 

i=l 

13: Delete Pk from list: list := list\Pfe; 

14: Append list by (P^, lb', ub'): 

If ti < glb_ub, then list := list U (P/., lb', ub'); If ub' < glb_ub, then glb_ub := ub' and P^ := P^; 
15: Append list by (P^, lb“, ub“): 

If t" < glb_ub, then list := list U (P^J*, lb“, ub“); If ub“ < glb_ub, then glb_ub := ub“ and 

_ TDU. 

16: Termination test: If glb_lb ~ glb_ub < e, then terminate := 1; Output P^ and its solution. 

17: Choose (P, lb,ub) from list such that the the associated lower bound lb is the smallest 

in the list and set the global lower bound glb_lb = lb. Let Pfe+i := P. 

18: k k P Ij 

19: end while 


5 . Extensions to uncertain Vis and MPCCs. In this section, we consider two key generalizations 
of the uncertain monotone linear complementarity problem. In Section [4.3[ we extend this framework to the 
regime of affine variational inequality problems over polyhedral sets. Next, we demonstrate how our frame- 
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work can address a subclass of stochastic mathematical programs with equilibrium constraints (MPEGs) 
(cf. [25jb given by a stochastic quadratic program with (uncertain) linear complementarity constraints. 


5.1. Uncertain affine monotone polyhedral Vis. Two shortcomings immediately come to the fore 
when considering the model (3.3): 


(i) The set X is a cone; 

(ii) The underlying set is deterministic in that it is uncorrupted by uncertainty. 


In this subsection, we show that examining uncertain polyhedral sets can also be managed within the same 
framework. Specifically, we begin by considering an uncertain affine variational inequality problem over a 
polyhedral set of the form given by (1.5) wherein X(u) and F(x,u) are defined as 


X(u) = {x : C(u)x > b(u), a; > 0} and F(x, u) = M(u)x + qiu)^ (5.1) 


respectively. From m Prop. 1.3.4], X solves Yl{X{u), F(u)) if and only if there exists a vector A G 
such that 


0 < a; T M{u)x — Ciu^^X + q{u) > 0 
0 < A T C{u)x — b{u) > 0. 


(5.2) 


In short, when F{x^ u) is an affine map and X{u) is a polyhedral set, the affine variational inequality problem 
is equivalent to a linear complementarity problem over a larger space of primal and dual variables. This can 
be more compactly stated as the following monotone linear complementarity problem: 


0 < z 1. B{u)z + d{u) > 0, 


(5.3) 


where 


Vc(o) 0 ) [-bM)’ 

respectively. It is relatively easy to see that B{u) is a positive semidefinite matrix since z'^B{u)z = 
x'^M{u)x > 0 if M{u) is a positive semidefinite matrix. This allows for making the following tractabil- 
ity claim when lA =hl 2 - Naturally, we may also extend other statements drawn from the regime of uncertain 
linear complementarity problems but leave that for future work. 

Proposition 5.1. Consider an uncertain variational inequality problem denoted by { VI{X{u),F{»] u)y\ueu 
whereU = {u \ ||u ||2 < 1}, X{u) and F{x,u) are defined in (5.1), where 


M{u) = S'^{u)S{u),S{u) = ^uiSi + 5o,g = go + y^ 


UiQu 






C{u) = Co + ^ uiCub{u) = hf) + w € U. 






Then a robust solution of this problem is given by a solution to a tractable convex program. 
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Proof. Recalling that z = (x, A), the robust counterpart of (5.3) is given by the following: 


min t 


subject to z^ 

{B{u)z + d{u)) < t, 

Mu G U, 

(5.4) 

M{u)x — 

C{u)'^X + q{u) > 0, 

Mu e U, 

(5.5) 


C(u)x — b{u) > 0, 

Mu e u, 

(5.6) 


x, A > 0. 




We begin by considering constraint (5.4) which can be recast as as z^[B(u) + B(uf)z + 2z^d(u) < 2t. 


Consequently, this constraint can be reformulated as a constraint similar to (3.16). By applying Lemma 


3.5 constraint (5.4) can be shown to be a tractable convex constraint. Next, we proceed to show the 


tractability of constraint (5.5) which is equivalent to the following collection of n optimization problems 
where i = 1,..., n: 


L 

{[sf Sq + SQSl]^.x + [qi] - - [Cfli.X) Ui + ^[SfSm + S'^Si]i.xuiu^ + ^[Sf Si]i.xuj 


L/=i 


l<m 


/=1 


subject to ||u ||2 < 1. 


(5.7) 


Analogous to Theorem |3.7[ this problem can be rewritten as a linear matrix inequality and some linear 
inequlities. It follows that (5.5) can be rewritten as a collection of n linear matrix inequalities and a bunch 
of linear inequalitiess. Finally, constraint (5.6) can be rewritten as the following set of constraints: 


[Coli.a; - [6o]i + ^ w/ ([C/]*.a: - [6/]*) >0, Mu & U, 


i = 1, 


,n. 


1=1 


This set of semi-infinite constraints is equivalent to a finite set of convex constraints in the form of second 
order cone constraints, which is discussed in Example 1.3.3. from 0- a 

5.2. Uncertain mathematical programs with complementarity constraints. Over the last two 
decades, the mathematical program with equilibrium constraints (MPEGs) has found utility in modeling 
a range of problems, including Stackelberg equilibrium problems, structural design problems, bilevel pro¬ 
gramming problems, amongst others. A comprehensive description of the models, theory, and the associated 
algorithms may be found in the monograph by Luo et al. |25j . When the lower-level problem is given by 
a complementarity problem, then the MPEG reduces to a mathematical program with complementarity 
constraints (MPGG). We consider the uncertain counterpart of MPGG defined as follows: 


min f{x,y) 

subject to h{x,y)>0 (5.8) 

0 < y 1. F{x, y) > 0. 

The MPGG is an ill-posed nonconvex program in that it lacks an interior. In fact, standard constraint 
qualifications (such as LIGQ or MFGQ) fail to hold at any feasible point of such a problem. We define an 
uncertain MPGG as a collection of MPGGs given by 


{MPGG(/,/i,F)}„ei^, 
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in which /, h and F are parametrized by u where u GU: 

min f{x,y,u) 

subject to h{x,y,u)>0, (5.9) 

0 < 2 / _L F{x, y, u) > 0. 

We may then define a robust counterpart of this problem as follows: 


min max f{x,y,u) 

uGU 

subject to h{x,y,u) >0, 'iuGll (5.10) 

0 < 2 / _L F{x^ y, u) >0, Vu G U. 


This problem is a nonconvex semi-infinite program. By utilizing the framework developed earlier, we may re¬ 
formulate (5.10) as a finite-dimensional MPCC. Unfortunately, the semi-infinite complementarity constraint 
given by 


0 < 2/ -L F{x, y, u) >0, \/u Ghl 


need not admit a solution. Instead, we recast the uncertain complementarity constraint as the following: 

min t 

subject to f{x, y,u) <t gU 

h{x, y,u) >0 yu gU 


min max y'^F{x,y,u) 

V u^U 

y solves F{x, y, u) >0, yu GU 

y>o 


(5.11) 


A natural question is whether a low-dimensional counterpart of (5.11) is available. Under convexity assump¬ 


tions on /(x, y, u) and concavity assumptions on h{x, y, u) in x and y for every u, and some assumptions on 


the uncertainty set, tractable counterparts may be constructed for the first two constraints in (5.11). By the 


findings of the prior sections, under some conditions, a robust counterpart of an uncertain LCP can be cast 
as a single convex program. The following result presented without a proof provides a set of assumptions 
under which the lower-level problem can be recast as a convex program: 

Proposition 5.2. Suppose F{x,y,u) is an affine map given by F{x,y,u) = Ax -I- M{u)y q{u) and 


M{u) = Mq + J2i=i uiMi,q{u) = q, Ml 0, VZ = 0,..., L , then the third constraint of (5.11) can be replaced 
by the optimality conditions of a convex program ifU =Ui,U 2 orU^o- 


6 . Numerical results. In Section [6T| we compare the quality of the residual of non-robust solutions 
with their robust analogues on an example presented in Section The performance benefits of the presented 
branching scheme on a non-monotone problem are examined in Section [6.2| and we conclude with a case study 
on uncertain traffic equilibrium problems in Section |6.3| where we compare robust solutions with the ERM 
solutions investigated in the literature. 

6.1. Monotone uncertain LCPs. We consider the constructed uncertain LCP defined in Section[2]for 
which the solution is known a priori. Table [O] shows that the presented techniques allow for obtaining the 
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robust solution and this corresponds closely with the analytically available solution Furthermore, 

an arbitrarily chosen scenario-specific solution, given as leads to large deviation from the analytical 

optimal solution and significantly higher residual. 


n 


residual of 

ll^n-rob _ 

residual of 

10 

3.9e-08 

2.0e-07 

0.4e-|-03 

5.0e+07 

20 

4.7e-08 

3.6e-07 

0.7e-|-03 

l.Oe+09 

40 

1.8e-07 

2.2e-06 

1.6e+03 

4.3e-|-10 

80 

5.1e-07 

5.2e-06 

5.5e+03 

3.9+12 

160 

1.6e-05 

5.3e-04 

1.6e+04 

2.8e+14 


Table 6.1: Robust vs non-robust solutions 


6.2. Non-monotone nncertain LCPs. We now consider a non-monotone LCP{M(u), q{u)) whose 
robust counterpart is given by the following: 


min t 

subject to X^{uiSi — U 2 S 2 )X + {uiqi + U 2 q 2 )'^X 
{uiSi - U2S2)x + Uiqi + U292 

X 

where = (l,...,n), Si = e„e^ ^ 0, S '2 = 10"^ x B > Q, U = {{ui,U 2 ) \ 0 < Mi ,-«2 < 1}, i? is a 
randomly generated matrix with elements drawn from A/'(0,1), qi = —e„ and q 2 = ^{n^i) ^ S 2 en- It should 
be emphasized that our analysis allows for deriving the robust counterpart of this problem as a relatively low¬ 
dimensional nonconvex QCQP. In the absence of such an analysis, a direct approach would require solving 
an approximate nonconvex QCQP whose size is of the order of magnitude of the discretization. Table |6.2| 
provides a comparison of the performance of three solvers on the RC on a set of test problems for increasing 
matrix dimension n: (i) our branching scheme, (ii) the commercial global optimization solver baron; and (hi) 
the multi-start solver from Matlab. The results from Table [6?^ suggest the following. First, our branching 


<t, Vu € U, 
>0, Vu G W, 
> 0 . 


( 6 . 1 ) 


Size 

Branching scheme 


baron 

matlab 



n 

time(s) 

^branch 

Nodes 

Gap 

time(s) 

^baron 

time(s) 

^matlab 

ll^brancli ^barony 
l+llxtiaronii 

l+||i;baron|| 

6 

1.32 

0.0100 

35 

0.00 

0.61 

0.0098 

3.52 

2.54 

0.00 

0.33 

7 

2.21 

0.0035 

82 

0.00 

1.81 

0.0164 

4.54 

2.53 

0.00 

0.21 

8 

2.29 

0.1648 

87 

0.00 

0.87 

0.1682 

5.07 

1.38 

0.00 

0.33 

9 

14.19 

0.0072 

406 

0.00 

0.89 

0.0076 

6.40 

3.65 

0.00 

0.23 

10 

9.80 

0.0040 

254 

0.00 

1001.7 

0.0125 

6.93 

0.74 

0.00 

0.22 

11 

72.74 

0.0036 

893 

0.00 

2.94 

0.0026 

6.09 

14.31 

0.00 

0.22 

12 

83.73 

0.1998 

1539 

0.00 

1.79 

0.1990 

8.22 

10.69 

0.00 

0.13 


Table 6.2: Global optimization of nonconvex QCQPs: CPU 3.40Ghz RAM 16.0 GB 


scheme provides reasonably accurate solutions by comparing with the commercial solver baron, somtimes 
even better with respect to the optimal value z. Furthermore, the performance is significantly superior in 
terms of optimal value to the solutions provided by Matlab. Third, baron’s performance in terms of time 
is superior to that provided by our Matlab-based branching solver is not altogether surpising, given that it 
uses extensive pre-processing and has been developed on C/C-l—1-. 

6.3. Case study: Uncertain traffic equilibrinm problems. 
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LU2 


Fig. 6.1: 2-node traffic network 


2-node and 5 link network:. Consider the uncertain traffic equilibrium of the form described in Section]^ 
sourced from m- Suppose the associated network has two vertices A,B and five arcs Di, Z? 2 , £^ 3 , C^ 2 - 

Let ^ denote the flow over these five paths and T{u)^ 1 represent the travel associated travel times, where 

T(it) is an uncertain 5x5 matrix and t S K® is a constant vector. Suppose B represents the path-OD pair 
incidence matrix and d{u) G represents the uncertain demand. Let r represent the minimum travel time 
for each direction. Recall that the equilibrium point is given by a solution to the following: 


0 < x _L M{u)x q{u) > 0, 
where x = M(u) and q{u) are defined as 


yuGU 




.B = 


1110 0 ' 
,0 0 0 1 1 , 


= 


( 1000\ 
950 
3000 
1000 
Vl300/ 


and T(u) is defined as 



^ 40a(M) 

0 

0 

20/3{u) 

0 \ 


0 

60/3(m) 

0 

0 

2QI3{u) 

T{u) = 

0 

0 

80/?(u) 

0 

0 


8q;(u) 

0 

0 

80q!(m) 

0 


V 0 

4/3(u) 

0 

0 

100;3(u)/ 

^u(u — 

1 ) and /3( 

u) = u 

(2 - u). 

Suppose U 

= {^1,142 


,d{u) = 



260 - 100(a(u) 
170- 100(a(u) 


a(u) = 

windy, and a rainy day respectively. In an effort to compare the obtained solutions with that obtained from 
the ERM model we assume that these events occur with probability 5 , |, and I (Note that our model 
does not require a probability distribution). Corresponding to this problem, the ERM solution is denoted 
by X®™ while the robust solution is Furthermore, non-robust scenario-specific solutions are denoted by 


x^,x^ and x^. Table 6.3 compares the optimality and feasibility of such points with respect to the robust 
counterpart. In this table, the infeasibility function is dehned as max(e^ max(—M(u)x — q{u),0)) while the 

uGU 

complementarity residual is defined as: ma.'xx^{M{u)x q{u)). It is seen that the robust solution and x^ are 

uGU 

feasible for every u. Notably, the “sunny day” design is feasible but leads to a large complementarity residual. 
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solution 

infeasibility 

complementarity 


(0,260,0, 170,0,950, 1000) 

0 

4.251E+06 

x‘^ 

(159.2, 0.83, 0, 70, 0, 1000, 1000) 

250 

1.717E-I-06 

x^ 

(0,160, 0, 3.75, 66.25, 950, 1300) 

500 

2.228E+06 


(84, 84, 21,80, 20, 975,1000) 

166 

1.089E-I-06 

a:”*’ 

(117.7, 89.5, 52.8, 90.5, 79.5, 950, 1000) 

0 

1.S40E+06 


Table 6.3: Comparison across solutions 



U — Ui 

U — U2 

U — Us 

a;”*’ 

1.4e+5 

1 .8e+6 

1 .8e+6 


Inf 

6.5e+5 

Inf 


OD pair 

possible demand 

Xro 

XeRM 

Xl 

X2 

X3 

AB 

260,160 

260 

189 

260 

160 

160 

BA 

170,70 

170 

100 

170 

70 

70 


Table 6.4: Evaluation of G{x*,u) Table 6.5: Flow of each OD pair 


Note that if both the feasibility and the complementarity metric is zero, this implies that the solution is an 
equilibrium for every u € U. The robust solution minimizes the worst complementarity residual among all 
possible scenarios and from that standpoint, it is seen to be superior to the solution that minimizes the 
residual for the first scenario. Furthermore, x®™ might have a superior complementarity residual but such 
a solution may be rendered infeasible for certain realizations. Table compares the value of parametrized 
gap function G{x,u), dehned as 


G(x, u) = sup(x — y)'^ {M{u)x + q{u)). 

y>0 

The lowest value of G{x,u) is achieved by x®™ when u = U 2 - However, G(x® 
consequence of infeasibility. However, G{x^°'^,u) < oo for every u GU. In Table 
robust solution satishes demand requirements (ensuring feasibility) while the ERM solution may not satisfy 
demand for all realizations (leading to infeasibility). 

5-node and 7-link network:. We now consider a larger traffic network considered in Sectionj^with 7 links 


ui) = G(x®^®,M 3 ) = + 00 , a 
we consider how the 


6.3 


and 6 paths. Figure 6.2 represents a 7-link network with 6-paths sourced from m and A ^ D and A ^ E 
represent two origin-destination (OD) pairs. The OD pair H —>■ D is connected by paths pi = {1,3},P2 = 
{1,7,6},P3 = {2,6} while the OD pair A ^ E is connected by paths p 4 = {1,5},P5 = {l,7,4},p6 = {2,4}. 
The demand along every OD pair is denoted by d{u) G where u GlA while the link capacity is captured 
by the vector c{u) gBJ ,u GU. Let vector x S 77® denote the assignment of flows to all path from pi to pg 
and / G 77^ denote the assignment of flows to all links 1,..., 7. Then the relationship between x and / is 
presented by: / = Ax, where A = (^ij) is the link-path incidence matrix. The entry Sij is set at 1 if and 
only if link i lies in path j. Let B = {bij) denote the OD-path incidence matrix and bij = 1 if and only if 
path j connects the zth OD pair. In this case, the two matrices are given as follows: 

/ 1 1 0 1 1 0 \ 


A = 


0 0 1 0 0 1 

1 0 0 0 0 0 

0 0 0 0 1 1 

0 0 0 1 0 0 

0 110 0 0 

V 0 1 0 0 1 0 / 


and B = 


1110 0 0 
0 0 0 1 1 1 


( 6 . 2 ) 
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Fig. 6.2: Traffic Network 


Following a generalized bureau of public roads (GBPR) function, the multivalue link cost function C{f,u) 
is defined as: 


a{f,u)=c^, 1.0 + 0.15 


i{u) 




(6.3) 


where c® and rii are known parameters. Let rii = 1 for all i, then the travel cost function is given as: 


F{u, x) = r]A^C{Ax, u) = O.lSr^A^diag ( ^ Ax + rjA^c^ A M{u)x + q. (6.4) 

\Ci{u)J 


Let w e denote the minimum travel cost of each OD pair. Last, by Wardrop’s user equilibrium, the 
uncertain CP formulation is given by the following: 


0 < 





Vu € U. 


(6.5) 


Tables 6.6 and 6.7 show the comparison between different solutions of the LCP given by (6.5). x^°^ denotes 
the robust solution of (6.5) in that it minimizes the worst case of the gap function G(z,m), defined as 


G{z, u) = sup( 2 ; — yy {A{u)z + b{u)),z = 

v>0 



,A{u) 


/m(u) 

V ^ 




We consider a case when U = {u \ —1 < u < 1}, —= (co)i +u(ci)i, where Cq and Ci are defined as follows: 

Cl \^U) 

Co = -Cl = (1/40,1/40,1/20,1/20,1/20,1/20,1/20), c° = (3, 5,6,4,6,4,1), 
d{u) = do + udi,do = (200; 220), di = (50; 40). 


The ERM solution xerm is constructed as follows. Let xerm = (y; w) where y is obtained by 

y = {I — B)x* + i?iE[d(M)], where B^ = B"’"{BB^)~^ . 

Note that x* is a minimizer of (f>{x) over the set H, where 

cj){x) = E[/(a;, u)], f{x, u) = z{x, uyF{z{x, u),u) + Q{z{x, u), u), z(x, u) = (I — B^B)x + B^d{u), 

F{z, u) = M{u)z + q, Q{z, u) = max ?/^E(z, u) \ By = d{u), y > 0} = min {^y'^d{u) \ B^y + F{z, u) > 0} , 

D = < X \ B"^Bx < c,Ci = min(i?id(u))i i , 

I u&A I 
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as per the recent work by Chen, Wets and Zhang m- Note that an estimator the minimizer of E[/(x,u)] 
is obtained via sample-average approximation schemes while w is acquired by taking the minimum of the 
average costs of paths in each OD-pair. Let the average costs of paths be captured by a vector v = 
¥.[M{u)y + q], then wi = min{ui, U 2 , U 3 }, 102 = inin{u 4 , U 5 , ue}. Xi, ■ ■■ are the solutions of the program 
iiiiiiA(ii)a;-i-b(«)>o,a:>o{®^(^('^)2^ + ^('^))}: when u = —1, —0.5,0,0.5,1, respectively. Table ( 


6.6 


shows the traffic 


OD pair 

range of possible demand 



Xi 

X2 

X3 

X4 

X5 

AD 

150-250 

250 

200 

150 

175 

225 

250 

200 

AE 

180-260 

260 

220 

180 

200 

240 

260 

220 


Table 6.6: Flow across two OD pairs 


flow between two OD pairs. Again, the robust solution satishes the largest possible demand while the ERM 
solution does not satisfy demand for all possible realizations. When we compare the residual function for a 
particular x and u as seen in Table |6.7[ while the robust solution x^°^ does not provide the best function 
value for every scenario, it minimizes the worst case. In fact, for the non-robust solutions, except x'^, every 
solution displays an infinite residual function for some u. Notably, the ERM solution also have infinite 
residuals for some realizations of u. 


u 

x"” 


xi 

x" 

x3 

x^ 


-1 

10343 

4340 

6488 

7922 

18322 

19000 

4329 

-0.5 

7863 

2176 

Inf 

6772 

14671 

14250 

2165 

0.0 

5382 

Inf 

Inf 

Inf 

11021 

9500 

0.000449 

0.5 

2901 

Inf 

Inf 

Inf 

7370 

4750 

Inf 

1.0 

421 

Inf 

Inf 

Inf 

Inf 

1.170e-05 

Inf 


Table 6.7: Residual function value at different sample points 


U 

M, q{u) 

M{u), q 

M{u), q{u) 

A'^ (u)A{u), q{u) 


convex QP 

convex QP 

nonconvex QCQP 

/ 

Hi 

convex QP 

convex QCQP 

nonconvex QCQP 

/ 

U-2 

convex QCQP 

convex program 

nonconvex program 

SDP 

Wc 

QP with conic constraints 

/ 

/ 

/ 

Ul 

/ 

convex QCQP 

/ 

/ 


/ 

convex QP 

/ 

/ 


Table 7.1: Characterization of robust counterparts under varying assumptions. 


7. Concluding remarks. In this paper, we consider the resolution of finite-dimensional monotone 
complementarity problems corrupted by uncertainty. A distinct thread in the literature has considered the 
minimization of the expected residual. This avenue relies on the availability of a probability distribution 
and the solution of a stochastic, and possibly nonconvex, program. Instead, we consider an avenue that 
relies on the availability of an uncertainty set. By leveraging findings from robust convex programming, we 
show that uncertain monotone linear complementarity problems can be tractably resolved as a single convex 
program. In fact, when the uncertain linear complementarity problem is not necessarily monotone, under 
some conditions on the uncertainty set, the tractable robust counterpart of this problem can be shown to 
be convex, a consequence of leveraging the hidden convexity in the problem. More generally, the robust 
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counterpart is a nonconvex quadratically constrained quadratic program. We adapt and present a recently 
presented branching scheme to accommodate such problems. Table |7.1| provides a compact representation 
of the tractability statements and the nature of the uncertainty sets that correspond to these statements. 
The columns of this table correspond to different assumptions of uncertainty on M and q. Note that 
M{u) = Mo + Y^f^^uiMi,Mi h 0 or ^ 0 for Z = 1,..., L, q{u) = qo + J2i=i'>J-iQh and A(m) = Ao + Ylf^^uiAi. 
We further observe that such statements can be utilized to show the following: 

(1) The tractable robust counterparts of an uncertain affine variational inequality problem(uncertain 
AVIs) over uncertain polyhedral sets are SDPs under some assumptions on the uncertainty set. 

(2) Robust counterparts of mathematical programs with uncertain linear complementarity constraints 
(uncertain MPCC) can be reformulated as deterministic low-dimensional mathematical programs 
with complementarity constraints. 

Future research will consider a study of nonlinear generalizations of F(x,u) as well as extensions to 
variational and hierarchical regime. Finally, our preliminary numerical investigations reinforce the belief that 
non-robust solution may produce large worst-case residual compared with robust solutions. Furthermore, we 
present a branching-based procedure for obtaining global solutions to robust counterparts of non-monotone 
uncertain LCPs and note its effectiveness compared to commercial global solvers. Finally, robust solutions 
are qualitatively different from their ERM counterparts in the context of traffic equilibrium problems. 

Finally, we comment on the nature of the uncertainty sets. We employ an “uncertain-but-bounded” 
model of uncertainty (a terminology that has its roots in in which the values of the uncertain parameter 
are represented through an affine parametrization of ^ which varies in the perturbation set Z. Naturally, this 
is by no means the only way to represent uncertainty. For instance, one alternate approach is to introduce 
chance constraints with ambiguity; in this model, we impose chance or probabilistic constraints under the 
caveat that the distribution is known partially in that it belongs to a family of distributions. We leave such 
questions for future work. 
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